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arbitrary linear combinations of the two modes in the definition of 
quadrature variables, we find that there always exists a combination 
of the two output modes which exhibits quadrature squeezing with 
noise reduction below the vacuum level. The number of noise photons 
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number of Stokes noise photons. 
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1 Introduction 



In this paper, we investigate novel aspects of quantum noise reduction in a 
model of a linear multi-channel quantum amplifier. Our investigation deals 
with both state reduction, or decoherence, and squeezing for the amplified 
channels. In fact, both phenomena turn out to be intimately connected. 
The model describes stimulated Raman scattering (SRS) in a loss-less cavity 
with a single Stokes mode, a single anti-Stokes mode and a classical pump. 
We assume the Raman medium to be not saturated, but allow for arbitrary 
homogeneous broadening and dispersion. The Raman amplification process 
gives rise to both two-mode squeezing and decoherence for the two optical 
field modes. We find that quantum noise resulting from these processes can 
be largely cancelled by appropriate interference arrangements, while main- 
taining the full amount of amplification. 

The first arrangement is an 'anti-squeezing' operation, by which two un- 
correlated normal modes are obtained. This operation can be implemented, 
for example, by a linear parametric amplifier. If the initial state of fields, 
medium and reservoir is the vacuum, the state for one of the normal modes is 
the vacuum, while for the other it is a canonical ensemble of number states. 
The average number of "noise" photons for the latter mode is equal to the 
number of phonons in the medium and reservoir states. 

The second arrangement involves linear mixing of Stokes and anti-Stokes 
without additional squeezing. The state of the resulting mixed mode is a 
canonical ensemble of squeezed number states. For a suitable, optimal choice 
of mode coefficients, the variance of one of the quadrature phases [1] of 
the mode is below the vacuum level, in spite of the additional quantum 
noise resulting from decoherence of the mixed state. This result holds for 
all values of dispersion, line width and coupling constants. The variance 
of the conjugate quadrature is comparable to the variances of the Stokes 
quadratures. For large amplification, the number of squeezed photons in 
the mixed mode is proportional to the square root of the number of Stokes 
photons. This represents a significant reduction in the amount of decoherence 
for the mixed mode, as compared to the original modes. 

The use of general mode coefficients amounts to a generalization of the 
notion of two-mode squeezing. For conventional two-mode squeezing, both 
mode coefficients are chosen to be of equal magnitude. (This is the defini- 
tion reported by Loudon and Knight [1]. Caves and Schumaker rescale the 
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coefficients by frequency dependent factors, to correct for the frequency de- 
pendent relationship between photon and field operators [2].) In this case, 
the variances of both quadratures and the number of photons of the mixed 
mode grow in proportion to the number of Stokes photons, indicating a much 
higher level of quantum noise than for the optimal choice of mode coefficients. 

The reduction of quantum noise through squeezing has been studied ex- 
tensively [1,3]. A squeezed state of a harmonic oscillator is a state for which 
the variance in one of the conjugate variables is reduced at the expense of an 
increase in the variance of the other variable. (We will not discuss the more 
general notions of higher order squeezing which have been defined in terms 
of the variance of polynomials of conjugate variables [4,5]. In the context of 
quantum optics it has become customary to speak of conjugate variables as 
quadrature phases [1]). 

For example, the variance in the electric field of an optical mode may be 
reduced below its vacuum level at the expense of an increased variance in the 
magnetic field. This situation is desirable for high sensitivity interferometer 
experiments [1,6] and optical communications [1,7]. Such squeezed states of 
the electromagnetic field can be generated by non-linear optical processes, 
such as degenerate parametric amplification or four-wave mixing, in which 
photons are created in pairs. If the process is non-degenerate, i.e. if two 
photons at different frequencies are generated, squeezing is not found for 
each of the individual modes, but rather for linear combinations of the two. 
This is known as two-mode squeezing, and its experimental detection requires 
an additional coherent source at the sum- or difference frequency [1]. Two- 
mode squeezing is conventionally defined in terms of linear combinations with 
equal mode coefficients [1,2]. We will show in this paper that, for situations 
where decoherence occurs, this is not an optimal definition, and that it is 
useful to adjust the mode coefficients to compensate for different noise levels 
in the two modes. 

The generation of squeezed light involves an amplification process, and 
this is frequently associated with state reduction or decoherence. State re- 
duction (the transformation of a pure state into a statistical mixture) oc- 
curs due to the coupling of the observed modes to a number of unobserved 
or unobservable modes. The process of state reduction in linear quantum 
amplifiers has been studied for various types of models, and the following 
references represent just a small selection of papers which are of interest here 
[8,9,10]. Such models are of great interest for the description of dissipation in 
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quantum systems [9,10], and for the analysis of state reduction in quantum 
measurement theory [8,10]. Most studies refer to harmonic oscillator models 
involving a small number of observable variables, and a large or infinite num- 
ber of unobserved 'reservoir' oscillators. The time evolution of the amplifier 
can be described either by master equations for the observable degrees of 
freedom [10], or by exact solutions for the linear equations of motion of the 
complete system [8,9,10]. 

The exact solutions reveal that, for these models, the time evolution is 
a multi-mode squeezing transformation coupling the observed variables and 
the reservoir. The amount of squeezing in the observable variables depends 
on the type of model studied. If the coupling to the reservoir involves only 
one of the conjugate variables, an initial rapid buildup of decoherence and 
squeezing for that variable occurs on a very short time-scale, of the order 
of the reservoir coherence time [9,10]. In addition, decoherence for both 
conjugate variables increases on a much longer time scale, given by the dissi- 
pation rate of the reservoir. In quantum optics it is customary to employ the 
'rotating wave approximation' (RWA), in which the coupling involves both 
conjugate variables of the observable system in a symmetric fashion. In this 
approximation, single- mode squeezing is absent, and decoherence occurs on 
the longer, dissipative time scale [9,10]. The cavity SRS amplifier studied 
here is in this category. 

So far, studies of quantum amplifiers in the RWA have emphasized either 
the process of state reduction, or the generation of squeezed states. However, 
both aspects are intimately connected. The authors found in [11] that, for 
the simplest case of a single mode amplifier coupled to an infinite reservoir, 
the time evolution results in two-mode squeezing for the amplifier mode and 
a single, time-dependent reservoir mode. If only the amplifier mode is ob- 
served, no squeezing is found, and the density matrix contains a large number 
of pure components, proportional to the amplification factor. If, on the other 
hand, both the amplifier mode and the large-time limit of the special reser- 
voir mode are observed, the density matrix contains only a small number 
of squeezed pure components. Similar results were found for a model of a 
single mode laser, comprising a single radiation mode, an inverted medium 
and a reservoir. To the best of our knowledge, these connections between the 
choice of observed mode combinations, the degree of state reduction and the 
degree of squeezing have not been systematically explored so far. 

The linear cavity Raman amplifier studied here is of particular interest, 
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both because of the type of interactions found in this system, and as an 
object for experimental study. The medium is assumed to contain a large 
number of unsaturated Raman active molecules or atoms, and is modelled by 
a single oscillator mode. The coupling between Stokes and medium oscillator 
is of squeezing type. However, since the medium is not directly observable, 
two-mode variables constructed from these two modes have no experimen- 
tal significance. On the other hand, the coupling between anti-Stokes and 
medium conserves photon-phonon number, and the anti-Stokes mode can 
therefore be considered as a probe of the medium. In fact, the squeezing 
transformation occurring between Stokes and medium is reflected in corre- 
sponding correlations between Stokes and anti-Stokes. While the coupling 
between anti-Stokes and medium is in most cases stronger than the Stokes 
coupling, the influence of the anti-Stokes on the medium is reduced by disper- 
sion, which leads to a frequency detuning of the anti-Stokes. It is therefore 
of interest to study the behavior of the CRA over the whole range of relevant 
parameters, including dispersion, homogeneous line width and gain. Experi- 
mentally, all of these parameters can be controlled independently, by varying 
the pressure of the gaseous Raman medium, the type of medium and the 
pump intensity. 

Cavity SRS models with Stokes and anti-Stokes modes have been studied 
by several authors. Non-linear, micro-maser models of SRS in a cavity with a 
single three-level atom and without dispersion were studied by Law et al. [12] 
and Puri et al. [13]. The latter found noise reduction below the vacuum level 
with the conventional notion of two-mode squeezing for certain values of the 
number of passes and the passing time of the atom through the cavity. Linear 
oscillator models of the medium and reservoir had been studied earlier by 
Walls [14] and Perina [15] for arbitrary dispersion and for the limiting case of 
large line width, where the medium oscillators are eliminated adiabatically. 
These studies were mainly concerned with photon statistics, and give no 
explicit discussion of squeezing. The first detailed study of squeezing in this 
type of models was done by Chizhov et al., who studied quadrature and 
higher order squeezing for the case of zero dispersion in a cavity [16,17]. A 
free-propagation amplifier was studied by Yeong et. al. [18]. In the limit of 
large line width [16], quadrature squeezing was reported to be absent, while 
such squeezing was found for the case of small line width in [17] and [18]. 
These studies discuss conventional two-mode quadrature squeezing in terms 
of the variance of one of the two-mode quadratures. 
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Our present work studies the linear cavity SRS model for arbitrary values 
of line width and dispersion. (The case of zero dispersion is rather special, in 
that Stokes, anti-Stokes and the medium undergo damped oscillations. For 
non-zero dispersion, exponential amplification occurs.) Two-mode quadra- 
ture squeezing is defined in terms of the ratio of the variances of two con- 
jugate two- mode quadratures, with arbitrary mode coefficients. We give 
explicit expressions for the Wigner characteristic function of the entangled 
state of Stokes and anti-Stokes modes, and determine those values of the 
mode coefficients, for which maximal squeezing and minimal state-reduction 
occur. 

In section 2.1, we give the model equations for the amplitude operators, 
derive the general form of solutions, and calculate the variance of a general 
Hermitean combination of the two optical modes for the case where the initial 
state at time t=0 is the vacuum. As is well known, this uniquely determines 
the Wigner characteristic function, and thereby the Schrodinger state at later 
times t > 0. In section 2.2, we show that the reduced state of Stokes and 
anti-Stokes is a product of the vacuum and a canonical ensemble of number 
states for two uncorrected normal modes. These normal modes are obtained 
from the original optical modes by an 'anti-squeezing' transformation, and 
we show that this transformation can be realized by a linear parametric 
amplifier. In section 2.3, we study 'non-squeezed' linear combinations of 
the two optical modes, and generalize the corresponding notion of two-mode 
squeezing by admitting arbitrary mode coefficients. We show that, at any 
given time t, it is possible to choose optimal mode coefficients in such a way, 
that one of the quadratures has minimal variance below the vacuum level 
at that time. This choice also leads to minimal decoherence. In situations 
where the system shows exponential amplification, the number of photons in 
the optimal mode is proportional to the square root of the number of Stokes 
photons, indicating a much higher level of quantum coherence than for the 
Stokes mode. In section 2.4, we consider amplification of an initial coherent 
state of the Stokes mode. The coherent state is maximally amplified if it is at 
90 degrees relative to the minimal two-mode quadrature, and gives the same 
output as the vacuum if it is parallel to that quadrature. In the former case, 
amplitude fluctuations are enhanced, while phase fluctuations are below the 
vacuum level. 

In section 3.1 we give the general solution to the linear equations of motion 
for optical modes and medium in terms of the eigenvalues of the 3x3 matrix 
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of coefficients. We discuss the special cases of zero and large line width in 
sections 3.2 and 3.3, and the intermediate case of finite line width in section 
3.4. In particular, we consider the influence of dispersion on the long-time 
behavior of the system. For zero line width, there exists a critical value of 
dispersion, below which the system shows periodic or quasi-periodic behavior, 
and above which exponential amplification occurs. For non-zero line width, 
the system always shows exponential amplification at large times, except for 
zero dispersion, where the behavior is that of damped oscillations. 

In section 4, we summarize our findings and discuss some possible implica- 
tions for a deeper understanding of state reduction in more general systems. 

2 Model equations, state reduction and squeez- 
ing 

2.1 Model equations and characteristic functions 

Our model of a cavity Raman laser contains the following variables: a classi- 
cal pump at frequency a> , a Raman- active medium with transition frequency 
o> 3 , a reservoir of oscillators with frequencies u r , r = 4...N, modelling homo- 
geneous broadening of the Raman transition, a Stokes mode at frequency 
Ui = ujq — 0U3, and an anti-Stokes mode at frequency a; 2 = + ^3 + A. The 
parameter A is a frequency detuning of the anti-Stokes mode which results 
from dispersion in the Raman-active medium. Figure 1 shows a level diagram 
of the medium, including the Raman transition, two transient intermediate 
states and the three optical frequencies u> , u>i, and uj 2 . 

The medium is assumed to be only weakly excited, with most of its popu- 
lation in the ground state. It can therefore be modelled by a single harmonic 
oscillator. We write the Heisenberg operators b n (t), n = 1...N for all vari- 
ables involved as products of time-dependent phase factors and amplitude 
operators b n (t): 

b^t) = e-^-^hit) , 
~b 2 (t) = e-^ 0+ ^%{t) , 
h(t) = e-^Xit) (k = 3,...,N). 

The Heisenberg equations of motion for the amplitude operators in the 
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Markov approximation are [16,17]: 



b\ = K,b 3 , (1) 

b 2 = -iAb 2 - K 2 b 3 , (2) 

N 

b 3 = Kl b\ + K 2 b 2 -Tb 3 + F, F = J2 g s e- lAa %(0) , (3) 

s=4 

K = -g r e lAr % (r = 4..JV), (4) 



where k±, k 2 and g r are the coupling coefficients for the medium with the 
Stokes, anti-Stokes and reservoir variables. The frequency detunings are 
A = u 2 — u 3 — l)q and A r = uj r — cu 3 . The solutions of these equations take 
the form: 

N 

b\(t) = T u (t)a[ + J2T lk (t)a k , 

k=2 
N 

bi(t) = Tai^al + ^TMat (i = 2...N) , (5) 

k=2 

where a n = b n (0), n = 1...N, and the coefficients T nm {t) are c-number 
solutions of equations (1) through (4) which satisfy: 

T nm (0) = 1; n = m 
T nm (0) = 0; n^m 

Because equation (3) for b 3 involves only the initial values b r (0), r = 4. ..TV, of 
the reservoir variables, the coefficients Tik{t) for i, k — 1, 2, 3 can be calculated 
from equations (1), (2) and a modified version of equation (3), for which 
F = 0: 

b 3 = Kib\ + K 2 b 2 - Tb 3 . (6) 

In the following, we will study the coherence properties of the Stokes and 
anti-Stokes modes in the state \ip(t) > at time t. The state is uniquely 
characterized by its Wigner characteristic function depending on t and 
arbitrary complex coefficients a n : 

X (t,a n ,a* n ) = <m\e i[m+4,H ° )] m)> 

= <^(0)|e i ^ ( ' )+ ^ t( ' )] ^(0)> (7) 
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where the operator (f>(t) is defined by: 

N 

<f>{t) = a$(t) + 5>A-(f) 

3=2 
N 

= ft(t)a\ + Y,Pk(t)a k (8) 

k=2 

From equation (5), the time-dependent coefficients (3 n are given by: 

N 



Pt(t) = a* 1 T 11 (t) + y £ t a i T il (t) 

i=2 
N 

k (t) = alT lk (t) + J2^iT ik (t) (k = 2...N) (9) 

i=2 

From the canonical commutation relations for the amplitude operators, one 
finds: 

N N 

muHt)] = -Iftl 2 + E l&l 2 = -l«il 2 + E N 2 (io) 

k=2 j=2 

If the initial state |-0(O)> is a product of vacuum states for all amplitude 
operators a n , the characteristic function x takes the form 

X (t,a n ,a* n ) = e~^ 

where 



<p 2 (t) = <m\m+^(t)) 2 \m> (n) 

V V 



El/3n| 2 = 2|A| 2 -| ai | 2 + E 
n=l i=2 



|2 



Here we have used equation (10) to eliminate the dependence on all coeffi- 
cients (3i with i > 1. If only the Stokes and anti-Stokes modes are observed, 
we set ojj = for i > 2, and obtain 

^ 2 (t) = 2|A| 2 -| ai | 2 + |« 2 | 2 (12) 
= ^Tn(t)+a 2 T 21 (t) 
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Equations (11) and (12) for the characteristic function of the two optical 
modes are the main results of this subsection. They show that the reduced 
state of these modes is completely determined by the two complex coefficients 
Tn and T 2 i. For observation of the Stokes mode alone, one sets a 2 = and 
finds 

V 2 (t) = (2|Tn| 2 -l)K| 2 = (2ni + l)|ai| 2 (13) 
For observation of the anti-Stokes mode alone, one sets a.\ — and finds 

cp 2 (t) = (2|T 2 i| 2 + l)|« 2 | 2 = (2n 2 + l)\a 2 \ 2 (14) 
Here is the average photon number: 

Ul = <m\bik\m> = { |T |^[| 2 1 

The characteristic function (11) of each of these modes is a canonical ensem- 
ble of number states, which is uniquely characterized by rii (see Appendix 
A.l). From equation (10), the following relationship is derived for the coef- 
ficients Tji 

N 

|Tn| 2 -|T 2 i| 2 = £|r fcl (*)| 2 = i2 2 >l (15) 
In terms of photon numbers, this is equivalent to 

N 

m-n 2 = J2 n k = r2 -1>0 (16) 

Ai=3 

showing that the difference between Stokes and anti-Stokes photon numbers 
is equal to the total number of medium and reservoir phonons. 

2.2 "Anti-squeezing" , Normal Modes and Decoherence 

For t > 0, the Stokes and anti-Stokes modes are correlated. Equating the 
coefficients of aiQ; 2 in equations (11) and (12), we find: 

<#))| hb 2 |^(0) > = T*,T 21 . (17) 
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We now seek to write the two modes q as linear combinations of uncorrelated 
normal modes q: 

b\ = Uc\ + V*c 2 (18) 
b 2 = Vcl + U* Cl ; \U\ 2 - \V\ 2 = 1 

In terms of the normal modes q, the mode operator <f) of equation (8) takes 
the form: 

cj){t) = a \b\ + a 2 b 2 = 7i*4 + 72C2 (19) 
with complex coefficients 7$ given by 

7* = a* ± U + a 2 V, 

72 = a\V* + a 2 U* , (20) 

and 

I 1 2 1 I 1 2 I |2 1 I |2 

-|ai| + |a 2 | = -|7i| + I72I 
Equation (12) for the variance cp 2 suggests to choose 

U = R^Tu , V = R' 1 ^ . (21) 
With this choice, (p 2 is diagonalized: 

ip 2 = (2R 2 -l) 1 2 1+1 2 2 (22) 
and the characteristic function x(7i 5 12) of the two normal modes factorizes: 

X(7i,7 2 ) = e-t^M'e- 1 *™*- (23) 

This shows that the reduced state for the two normal modes is the product 
of a canonical ensemble of number states for c\ (see appendix A.l) and the 
vacuum state for c 2 . In particular, the normal modes are uncorrelated and 
not squeezed. The number n[ of photons in mode c\ is equal to 

n[ = R 2 - 1 = m - n 2 . (24) 

For the special situation where R 2 — 1, the coefficients Tn for i > 2 vanish. 
The normal modes q are then just the Schrodinger operators q(0) = c^, and 
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equation (18) is the original time-evolution, equation (5). As we shall see 
below, this situation occurs for T=A=0 periodically at certain times t > 0. 

It is remarkable that, even for R 2 > 1, one of the normal modes is always 
found in the vacuum state. 

A physical interpretation of the normal modes q can be given in terms 
of an "anti-squeezing" experiment, performed on the output of the cavity 
Raman laser. Inverting equations (18), we find 

c\ = U*b\-V*b 2 

c 2 = -Vb[ + Ub 2 (25) 

This is a squeezing operation performed on the modes b\ and b 2 , which can be 
implemented, for example, by a non-degenerate linear parametric amplifier 
(LPA) [3]. If the Stokes and anti-Stokes modes are used as idler and signal 
in a phase-matched linear amplifier, the equations for propagation in the 
z-direction are [3]: 

ge~^c 2 ; c\(0) = b\ 

ge^c\ ; c 2 (0) = b 2 (26) 

where 77 is the pump phase, and g is the product of pump amplitude and 
gain coefficient. The solutions are: 

c\(z) = cosh gz b\ + e~ ir) sinh gz b 2 

c 2 (z) = e tri sinh gz b\ + cosh gz b 2 (27) 

Equation (25) is obtained, apart from a common phase factor of c\ and c 2 , 
for 

cosh gz — \U\ ; rj — arg(V) — axg(U) (28) 

Equations (19) and (23) show that, at any given time t, this operation will 
disentangle the state of Stokes and anti-Stokes modes, leading to zero corre- 
lation between the new modes C\ and c 2 and to zero photon count for mode 
c 2 . We illustrate the sequence of squeezing and "anti-squeezing" operations 
performed by the cavity Raman laser (CRL) and the LPA in figure 2. 
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By adjusting the pump amplitude and phase of the LPA, arbitrary "anti- 
squeezing" transformations can be performed on the optical fields. For ex- 
ample, if the pump phase r) is changed by an amount 2e from the value given 
in (28), we find for the number of photons in the exit mode c 2 (z): 

T] = arg(F) - arg(C/) + 2e (29) 

and 

< V(0)| cl(z)c 2 (z) |^(0) > = 2 ( ni + 1 ^ 2 sin 2 e >2n2 sin 2 e. (30) 

1 + ri\ — n 2 

For n 2 ^> 1, a single photon is counted for e=(2n 2 )~^, and a very strin- 
gent condition on rj is obtained. Similar restrictions are found for the pump 
amplitude g. 



2.3 Generalized two- mode Squeezing for Stokes and 
anti-Stokes 

Two-mode squeezing for two oscillators b\ and b 2 is conventionally defined in 
terms of the variance of the Hermitian components [1]: 

X Q = ±(b 1 + b\ + b 2 + bl) 

X P = y(6i-6l + 6 2 -4) (31) 

These operators are just the quadratures of a linear combination b of bi and 
b 2 : 

b = ^fc + ^fe 02) 

X Q = ^(& + & f ); X P = ^(6-6t). (33) 

A state is two-mode squeezed for b\ and b 2 if it is squeezed for b. The con- 
ventional definition of squeezing is that the variance of one of the Hermitian 
components is below the vacuum value \: 

< (AXp) 2 > < \ for F = P or Q . (34) 
13 



Because of the uncertainty relations, this implies that 

<(AX P ) 2 >^<(AX Q ) 2 > . (35) 

In this paper, we use the term "squeezed state" for states satisfying the 
weaker condition (35). If the state is a minimal uncertainty state, e.g. the 
vacuum state, one has 

<(AX P ) 2 x(AX Q ) 2 > = ^ (36) 

and equation (35) implies (34). For states satisfying equation (35) but not 
(36), however, the variance of both components may exceed the vacuum 
value. 

To generalize the notion of two-mode squeezing, we first admit arbitrary 
coefficients Zj in the definition of the linear combination b 

b = z 1 b l + z 2 b 2 ; |zi| 2 + \z 2 \ 2 = 1 . (37) 

Second, we allow arbitrary phase factors in the definition of quadratures. We 
thus consider a general Hermitian linear combination of b and as in the 
definition (7) and (8) of the characteristic function 

fa = a*b ] + ab (38) 

For a = ^ or a = one has 4>i — Xq or fa = X P . To determine the 
degree of squeezing of a given state, one needs to find those values of a with 
\a\ = 4g, for which the variance of fa assumes maximal and minimal values 
ip 2 + and tp 2 _ . We now apply this procedure to the Stokes and anti-Stokes 
modes defined above. 

From equations (38) and (8) we find that the characteristic function of 
the mode b in equation (37) is obtained from the general expressions (11) 
and (12) by setting a\ = az\ and a 2 = az 2 . Thus 

(3{ = c?z\T xl + az 2 T 21 . (39) 

This gives for the variance kp\ of 4>\ (see equation 12): 

iff = 2\a*z*T n + az 2 T 21 \ 2 + H 2 (|z 2 | 2 - N 2 ) . (40) 
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To quantify the amounts of state reduction and squeezing, we rewrite this in 
the form (see Appendix A.l) 

<pl = Q 2 \ a *U' + aV'\ 2 ; \U'\ = cosh s , \V \ = sinh s. (41) 

This corresponds to a canonical ensemble of squeezed number states with 
average photon number n = \{Q 2 — 1) and squeezing parameter s. For 
| a | = we find the maximal and minimal variances of the corresponding 
quadratures from equations (40) and (41): 

*4 = (nRi ± r 2 R 2 f + \{r\ - r\) = \q 2 e ±2s , (42) 



where r\ = \zi\ and R% = \Tn\ . (43) 
This gives for the coefficients Q 2 and e~ 2s 



Q 2 = 2^pLrt , e~ 2s = yW; 2 (44) 
We now consider a situation where 

R 2 

Ri,R 2 — > oo , — > q < 1 as t — > oo . (45) 

Ri 



and give the results for various choices of mode coefficients Zj. 
C 

gives 



Case 1: The conventional choice of mode coefficients is \zi\ = ^ and 



Vl= l -{R 1 ±R 2 ) 2 ^ 1 -R 2 1 {l±q), 



Q 2 = R 2 1 -R 2 2 ^R 2 1 (l-q 2 ), ^^YVq (46) 

In this case, the quantum noise in both quadratures increases as R\ , and goes 
above the vacuum level. The average squeezed photon number increases in 
proportion to R\ , while the squeezing coefficient goes to a finite limit. 

Case 2: If the mode coefficients Zj are adjusted at any given time to make 
the first term in the expression for ip 2 _ in equation (42) vanish, we obtain: 

\zlTnl = \z 2 T 2l \ , (47) 
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2 1/2 2\ 1/1 9 \ . 1 2 i 2n2 a r>2 Q 

f- 2 v 2 i; 2 v l + g 2; 2 ' ^ + 11 1 1 + ? 2 

In this case, the quantum noise in the quadrature corresponding to y? 2 stays 
below the vacuum limit. The average squeezed photon number increases in 
proportion to Ri, while the squeezing coefficient decreases in proportion to 
R~,\ 

Case 3: It is possible to choose the mode coefficients Zj at any given time 
such as to minimize the variance (p 2 _. This minimal value is given by: 

vLn = \{r\+ri- -ri- i) 2 + < 2 («) 

In the limit of large times where equation (45) holds, this gives the same 
results as for the choice of coefficients in equation (47) above. 

Case 4: It turns out that, in many situations, noise reduction below the 
vacuum level can be achieved with constant mode coefficients and quadra- 
tures. For the choice a — we get for the variance in this case 

<p% x = |z?r u - z 2 T 21 \ 2 + \{\z 2 \ 2 - N 2 ) . (49) 
2.4 Amplification of coherent states 

In this section, we calculate the final state for the linear combination mode, 
as discussed in section (2.3), for the case where the initial state is a coherent 
state for the Stokes field, and the vacuum state for all other variables. 

ai|V(0)> = vi |^(0) > , ai\ip(0)> = for i > 1 (50) 

In this case, the characteristic function takes the form 

X (t,a,a*) = e i ^ Vl+ ^- 1 ^ . (51) 

If we take the near-optimal choice of case 2 for the mode coefficients z i: for 
which 

I^Tnl = \z 2 T 21 \ = ni?i , (52) 
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we find for the phase of the characteristic function (51) from equation (39) 

to + to* = (aS + a* 6*) with 6 = (i^ + vDnRx = 6* . (53) 

For the variance <p\, we obtain from equation (40) 

ip\ = 2\a* + a\ 2 r 1 R l + \a\ 2 {r 2 2 -r 2 1 ) (54) 

The characteristic function is that of a canonical ensemble of coherent, squeezed 
number states with coherent displacement parameter 6 (see Appendix A.l). 
The parameter 8 is maximal for real ui, and vanishes for imaginary v\. On the 
other hand, the variance f \ is maximal for real a, and minimal for imaginary 
a. The amplified state therefore shows phase squeezing. The amplification 
and squeezing transformations in the complex v plane are illustrated in figure 
3. 

3 Solutions of the model equations 

3.1 General case 

The matrix coefficients (t) are obtained as the classical solutions to equa- 
tions (1), (2), and (6) with T n (0) = 8 n . They take the form: 

= 5> ifc e A *' (i = 1,2,3) (55) 

k 

where are the eigenvalues of the matrix representing the right-hand side of 
equations (1), (2), and (6) (see appendix A. 2). The general form of the eigen- 
values depends on the relative size of the gain coefficients «j (i — 1, 2). These 
coefficients contain a factor involving the frequency of the mode considered. 
Except for a situation where the Stokes coefficient shows large resonant en- 
hancement over the anti-Stokes coefficient, one has k 2 >k\. We restrict our 
discussion to this case. 

In this case, the eigenvalues are either purely imaginary, or there is an 
eigenvalue with positive real part. In the former case, the matrix coefficients 
Tn show periodic or quasi-periodic behavior in time. In the latter case, the 
term involving an eigenvalue with positive real part dominates (55) for large 
times, and equation (45) holds. We now discuss some special cases. 
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3.2 Hypertransient limit 

In the hypertransient limit, we set T = 0. There are two regimes, determined 
by a certain critical value A cr u of A. For A < A cr u, the anti-Stokes mode 
dominates. The eigenvalues are all purely imaginary and the system shows 
quasi-periodic behavior. For A = 0, one of the eigenvalues is zero, and the 
matrix coefficients Tn(t) are periodic, with frequency 




At times t such that cos nt = —1, we have 

|Tn| = i?i = k- 2 (k\ + k 2 2 ) , \T 21 \ = R 2 = 2k~ 2 k 1 k 2 , R\ - R\ = 1 

This is the special situation mentioned in section 2.2 above. In this case, 

minimal variance is found for conventional two-mode squeezing \zi\ = \z 2 \ = 
i . 

2 2 1 /«2 - ^l\ 2 /, 7 x 

fconv = fmvn = « ~ T ~ J I 57 ) 

In the following numerical examples we choose K\ — 1 and k 2 = v^2- Figure 
(4.1) shows c^ in and <£ 2 conv as functions of t for A = 0. 

Figure (4.2) illustrates the quasi-periodic behavior of ^> 2 min and (p 2 ix for 
A = 0.290, just below the critical value A crit = 0.300. In this case, two 
of the eigenvalues are very close to each other, (A3 — A2) = ie, |e| <C |Afc|, 
and the matrix coefficients Tn(t) are mostly dominated by the corresponding 
terms with small frequency denominator e. Consequently, the variance ^p 2 min 
is almost constant, except for short time periods, during which these terms 
are small due to their explicit time-dependence on sinet. The behavior of 
<Pf ix is quite erratic in this case. The variance tp 2 conv is much greater than |, 
except for the short time periods mentioned, and is not shown in this figure. 

For A = A cr it, the matrix coefficients grow linearly in time, while for 
A > A crit , one of the eigenvalues has a positive real part. As mentioned 
above, the corresponding term in equation (55) dominates the behavior for 
large times, and the variance <p 2 min approaches a constant value. For this 
situation, we also plot the variance p 2 ix obtained for a fixed choice of mode 
coefficients Zi, which gives almost minimal variance for large t. Figure (4.3) 
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shows tpfnini f}ix an d ^lonv f° r A = 0-4, just above the critical value, while 
figure (4.4) shows the same for A = 10, much larger than A crit . In the former 
case, (fl onv is close to <^ 2 min and less than | up to t — 4, while in the latter 
case fl onv > | at all times. In both cases, v?„j„ and (p 2 ix go to different 
values (less than |) as t — > oo. 



3.3 Steady-state case 

In the steady-state limit, it is assumed that the homogeneous line width V 
is much larger than any of the frequencies K\, k 2 and A. In this limit, the 
mode operator b 3 is approximated by its steady-state value, and we obtain 
the following equations: 

b 3 = r- 1 (K 1 b\ + K 2 b 2 + F) 

b\ = 9iib{ + + Ci (58) 
h = -gi 2 b\ - g 22 b 2 + «A6 2 + G 2 

where 

Gi = KiY^F and gij = T~ 1 K i Kj 

These equations are equivalent to the model studied in references [14] and 
[16]. The matrix coefficients Tn are given by a formula analogous to (55), 
but involving only the two eigenvalues corresponding to equation (58). The 
eigenvalues are: 



A = ~(<?22 - gu - iA) ± ^(<722 - 9ii) 2 - ^A 2 - l -iA{g 22 + g u ) (59) 
For A = 0, one has: 

Ai = 0, A 2 = -(g 22 - g n ) (60) 

and the matrix coefficients have finite limits for t — > oo. In this case, one 
finds: 

Rl — > K~ 2 k\ , R 2 — > K~ 2 K\K 2 
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Figure (5.1) illustrates this case. 

For A > 0, one of the eigenvalues is positive , while the other is negative. 
In this case, the matrix coefficients diverge for t — > oo, and ^ 2 min goes to the 
finite limit given by equation (47). Figures (5.2) and (5.3) show the cases 
A = 0.1 and A = 1. They demonstrate that y> 2 min and ip^ ix converge to the 
same limit. The conventional variance fl onv exceeds the vacuum value \ for 
sufficiently large time. For A ^> g^, the eigenvalues are given by: 

Ai = gu and A 2 = -#22 + *A (62) 

and the Stokes and anti-Stokes decouple in this limit [14]. 

3.4 Transient case 

In the transient case, where T is of the same order of magnitude as K\ and 
re 2 , the cases A = and A>T are easily solved. For A = 0, the eigenvalues 
are: 

r Ft 2 

A = 0, \ = --±^--tf (63) 

and the matrix coefficients converge to the same limits as in the steady-state 
case. For small T, two eigenvalues have imaginary parts, and additional 
oscillations occur for small times. For A ^> T, the eigenvalues are, to leading 
order: 

A = *A, \ = - T -±^ + k\ (64) 

One eigenvalue has positive real part, and the mode coefficients grow expo- 
nentially for large times. For V ^> k\ , the results of the steady-state limit are 
recovered. For intermediate values of A, one eigenvalue has positive real part, 
and the behavior of the variance for large times is similar to the steady-state 
case. 

4 Summary and Discussion 

For the model of a linear cavity Raman amplifier studied here, we found 
two main results. First, there exist uncorrelated normal modes, which are 



20 



obtained from the Stokes and anti-Stokes modes by an anti-squeezing trans- 
formation. Second, there exist optimal linear combination modes, for which 
the variance of one quadrature is always below the vacuum value. The states 
of such optimal modes are canonical ensembles of coherent squeezed number 
states, and the average number of noise photons is proportional to the square 
root of the number of Stokes noise photons. Thus, the degree of decoherence 
for the combination mode is much less than for the individual Stokes and 
anti-Stokes modes. 

The reduction of variance in the minimal quadrature of the optimal com- 
bination mode can be understood as resulting from a cancellation of anti- 
correlated quantum noise in the individual modes. To optimize the amount 
of cancellation, the mode coefficients have to be chosen appropriately, and 
will in general not be of equal magnitude. This corresponds to a generaliza- 
tion of the conventional notion of two-mode squeezing. 

If a coherent state is used as input to the Stokes channel of the amplifier, 
the corresponding field amplitude will be cancelled in the combination mode, 
if it has the same phase as the minimal quadrature, and will be optimally 
amplified if its phase differs by ±90 degrees. The optimally amplified signal 
is thus phase-squeezed. 

These results may be of interest for the generation of squeezed light 
through stimulated Raman scattering. They also provide clear experimen- 
tal signatures for the experimental observation of anti-correlation between 
Stokes and anti-Stokes. It is of particular interest, that these signatures are 
found for all values of the experimentally relevant parameters, such as pump 
amplitude, dispersion and line widths. Previous work on quadrature squeez- 
ing in this model had been restricted to a limited range of parameters (small 
linewidth, zero dispersion). 

With the recent demonstration of a cw-Raman laser by J.K. Brasseur et 
al. [19], experiments to verify our predictions may become feasible in the 
near future. 
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Appendix A.l 

The quantum-mechanical density matrix p of a system with creation and 
annihilation operators at and a is uniquely characterized by its Wigner char- 
acteristic function 

X (a, a*) = Trace (pe^ (a ' a * } ) = < e^ (a ' a * } > (65) 

where 

4>(a, a*) = aa + a*a) . 
Of special interest are Gaussian states, for which 

x = e i«t»-\<{^?> ) A0 = -<0> . (66) 

The most general Gaussian state is a canonical ensemble of coherent squeezed 
number states. Using the notation of Caves [1], we write this as 

p = D(6)T(C)p th T- 1 (OD- 1 (6) 
with D(6) = e Sa ~ 6 * at , T(C) = e^ a2 ~^ an , ( = se ld (67) 
and p th = (1 - e~ r ) ^ e~ rn \n >< n\ . 

n 

For the expectation value and variance of <f>, one finds [1]: 

<</)> = ad + a* 6* 
<(A(f)) 2 > = Q 2 \aU + a*V\ 2 (68) 

with: 

,2 „ < l + e- r t e~ r 

) — 9n J- 1 — ■ n — /it 



Q 2 = 2n + 1 = - — ; n=<a T a> = 



1 - e~ r 1 - e" r 

[/ = cosh s, V = -e w sinh s (69) 

Equations (66), (68) and (69) give the characteristic function in terms of the 
coherent-state displacement parameter 6, the squeezing parameter ( through 
U and V and the mean photon number n through Q as: 

X (a, a*) = e iW+a*6*)-hQ 2 \«u + a*v\* (7Q) 
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Appendix A. 2 



The matrix coefficients Tn(t) for % — 1,2,3 satisfy equations (1),(2) and (6), 
with initial conditions 



Tn(0) = 1 and T a (0) = for i = 2,3 
They take the form 

^i(t)=E^ eAfc * * = 1,2,3 



(71) 



(72) 



where A& are the eigenvalues of the matrix of coefficients iV for equations 
(1),(2) and (6): 

/ «! 
7V = zA -«2 | (73) 



We find: 



(XjXj + k\) 



T2k 



d k 



i4(T-S k ) 
d k 



(74) 



where 



dk = (A fc - Ai)(Ajfc - Xj) and S k = Xi + Xj 

and, for fixed k, the indices % and j take the two values different from k. The 
eigenvalues were obtained numerically. 

Analytical solutions or approximations for the eigenvalues were obtained in 
the special cases considered in subsections (3.2) and (3.3), in particular the 
cases: T = 0, r 3> A, nf , A = 0, A 3> T, k? , and were used to verify the 
correctness of the numerical procedure. The calculations are straightforward 
and will not be detailed here. 
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Figure Captions 



Figure 1. Shown are the two levels of the Raman transition at frequen- 
cies uj = and uj = uj^, and two intermediate levels at frequencies uj = ujq 
and uj = ujq + 0J3 for Stokes and anti-Stokes generation. The corresponding 
absorption of pump radiation and generation of Stokes and anti-Stokes is 
indicated by vertical and horizontal arrows. The anti-Stokes mode at fre- 
quency uj = uj + uj 3 + A and its driving polarization at uj = uj + uj 3 are 
represented by two different lines, separated by the frequency detuning A. 

Figure 2 illustrates the sequence of entangling and dis-entangling operations 
performed by the cavity Raman laser (CRL) and linear parametric amplifier 
(LPA). The uncorrelated states of the initial Stokes and anti-Stokes modes ai 
and et2, and of the final normal modes c\ and C2, are represented by disjoint 
circles, while the entangled state of the intermediate, amplified Stokes and 
anti-Stokes modes b\ and 62 is represented by a single ellipse. The quantities 
rii and n 2 are the corresponding photon numbers. 

Figure 3 illustrates the effect of amplification on the coherent states of the 
near-optimal combination mode b. The coherent states are represented by 
ellipses in the complex v plane, which are centered at the coherent displace- 
ment parameter v = v\ for the initial states and v = 6 for the final states. 
The axes of the ellipses are proportional to the root of the variances of the 
corresponding quadratures. The arrows connect the initial state to the final 
state resulting from amplification. Case (a) represents the vacuum going to a 
squeezed state. Cases (b), (c) and (d) represent various initial coherent states 
with imaginary, real and complex 6 going the corresponding final states, with 
(b) going to the same state as (a), and (c) going to the same state as (d). 

Figures 4. a to 4.d. Shown are the variances of the minimal quadrature vari- 
able for the choice of conventional, fixed and minimizing mode coefficients 
(^conv; V'fix an d V^min) a ^ zer0 hnewidth (r = 0), for various amounts of dis- 
persive detuning as given in the figures (A = 0.0, A = 0.29, A = 0.4 and 
A = 10.0). 

Figures 5. a to 5.c. Shown are the variances of the minimal quadrature vari- 
able for the choice of conventional, fixed and minimizing mode coefficients 
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(v^conv) V'fix an d ^min) & t large linewidth (r = 5), for various amounts of 
dispersive detuning as given in the figures (A = 0.0, A = 0.1 and A = 1.0). 
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5 Figures 
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